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Abstract 

The problem is analyzed of extrapolating power series, derived for an asymptot- 
ically small variable, to the region of finite values of this variable. The considera- 
tion is based on the self-similar approximation theory. A new method is suggested 
for defining the odd self-similar factor approximants by employing an optimization 
procedure. The method is illustrated by several examples having the mathematical 
structure typical of the problems in statistical and chemical physics. It is shown 
that the suggested method provides a good accuracy even when the number of terms 
in the perturbative power series is small. 
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1 Introduction 



The problem of extrapolating power series is constantly met in many branches of physics 
and applied mathematics, when one considers realistic systems with strong interactions. 
The theoretical description of such systems can be reahzed by invoking some numerical 
techniques. For instance, the investigation of model fluids can be done by molecular dy- 
namics simulations, as is reviewed in Ref. [1]. However, it is always desirable to have 
analytical equations allowing for an explicit description of the matter properties. The 
standard way of obtaining analytical equations is by using some kind of perturbation 
theory. But the difficulty in deriving such analytical equations is that the realistic matter 
not often possesses small parameters, with respect to which one could develop perturba- 
tive expansions. Vice versa, the usual situation is the absence of such parameters. It is 
possible, nevertheless, to formally employ perturbation theory with respect to some sys- 
tem parameters, keeping in mind that the resulting series necessarily require an effective 
resummation. Among the methods, defining an effective limit of asymptotic series, the 
most popular are the Fade approximation [2] and the Borel summation [3]. But, as is well 
known, the method of Fade approximants is not uniquely defined and is plagued by the 
occurrence of annoying divergences. While the Borel summation is applicable only for the 
cases enjoying a large number of perturbative terms and the knowledge of the high-order 
expansion coefficients, which is rarely available. Recently, a novel general approach has 
been developed for extrapolating perturbative series, called the self-similar approximation 
theory [4-13]. It has been persuasively demonstrated by numerous examples [4-13] that 
this approach is much simpler and more accurate than other extrapolation techniques. 

The term self-similar with respect to the considered approximation theory is used 
because of the basic idea behind this approach. This idea can be described as follows. 
Suppose that by employing perturbation theory, we derive a sequence {fk} of approxi- 
mations fk, with the index k = 1,2... enumerating the approximation order. Considering 
the index k as the variable similar to discrete time, we can treat the motion from the 
A;-th approximation term to the higher-order terms. A; -|- 1, A; -|- 2, ... as an evolution of a 
dynamical system. A dynamical system with discrete time is called cascade. Then the 
sequence {fk} represents the points of a cascade, which form the cascade trajectory. If 
we would know the equation of motion for this dynamical system, then, starting from the 
initial points of the trajectory with low-orders k = 1, 2, we would be able to find the 
higher-order terms, including the point representing an effective limit /* of the sequence 
{fk} as k — > GO. This effective limit corresponds to the fixed point of the cascade tra- 
jectory. In the vicinity of the fixed point, it is possible to write down the equation of 
motion as a functional relation between the terms fk+p and fk, where p > 0. This relation 
has the form of the renormalization-group equation, which is also called the equation of 
functional self-similarity or the relation of functional self-similarity. Our approximation 
theory is based on the use of this self-similarity relation, because of which the method is 
named the self-similar approximation theory. All mathematical details and foundations 
of this approach can be found in Refs. [4-12]. 

A very convenient and powerful variant of the self-similar approximation theory is the 
method of self-similar factor approximants [14-18]. The direct use of this method requires 
the knowledge of even numbers of terms in the asymptotic series. However, it would be 
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certainly desirable to be able to work with any number of terms, whether even or odd, 
in order to use in full all the information contained in the given series. A simple way of 
avoiding this problem was employed in Refs. [16-18], where it was suggested to set one 
of the parameters of the odd factor approximants to unity, which could be explained by 
scaling arguments. 

In the present paper, we propose a more justified method for defining the odd factor 
approximants. This method is based on the optimization procedure, which follows the 
ideas of the optimized perturbation theory [19,20] (see also review articles [21,22]). We 
demonstrate by a number of examples that this new method provides a good accuracy 
of extrapolation for various problems typical of statistical and chemical physics. The 
method works well even when the number of the terms in the asymptotic series is rather 
small. Comparing the results of the calculations, employing the suggested optimization 
procedure, with the results of the method using the scaling arguments [16-18] makes it 
possible to justify the latter method. 



2 Formulation of problem 

Suppose we aim at finding a real function f{x) of the real variable x. But the equations, 
defining this function, allow us to derive its form only for asymptotically small x, where 
the sought function 

f{x):^Mx) (x^Q) (1) 

is approximated by the series 

fk{x) = fo{x)j2^nX''. (2) 
n=0 

Here /o(a;) is a term defined so that Gq = 1. 

The extrapolation of this asymptotic series to finite values of the variable x can be 
done by means of the self-similar factor approximants [14-18] prescribed by the expression 

Nk 

f;{x)^fo{x)l[{l + Axr, (3) 

i=l 

with the number of factors 

^ U-/2, /. = 2,4,... 

'''' \ (A; + l)/2, A; = 3,5,... " 

The parameters Ai and powers rii are defined by the accuracy-through-ordcr procedure 
that means the following. We expand the factor approximant (3) in powers of the variable 
X up to the k-th order. The resulting expansion is compared with the initial series (2) 
and the coefficients at like orders of x are equated with each other. Or we may compare 
the expansions of A;-th order obtained from lnfk{x) and ln/^(a;). Both these ways are 
evidently equivalent, since the logarithm of a function and the function itself arc in one- 
to-one correspondence with each other. The second way gives slightly more convenient 
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equations. Following this way, with equating the like-order terms in the expansions for 
In fk{x) and lnf^{x), we come to the system of equations 

^n,A^ = S„ (n = l,2,...,A;) , (5) 

1=1 

in which 

{n-l)\ 0.-^0 dx- ) 

All values of Bn are expressed through the coefficients ai, 02, a„ of expansion (2). 
Explicitly, Eqs. (5), e.g., for even k, are 

niAi + 712^2 + . . . + nk/2Ak/2 = Bi , 

niA\ + ^2^12 + ■ ■ ■ + "'fc/2^fc/2 = -62 , 



+ naA^ + . . . + nfc/2A^/2 = Bu . 

The solutions to these equations involve Vandermonde determinants. When k is even, 
there are k unknowns, k/2 parameters Ai and k/2 powers n^. Then all these unknowns 
are uniquely defined. But when k is odd, there are k -\- 1 unknowns, that is, (A; + l)/2 
parameters A^ and {k + l)/2 powers n^. Then, to define uniquely all parameters, it is 
necessary to impose an additional condition defining, for instance, Ai. In Rcfs. [16-18], 
we used the condition Ai = 1, following from scaling arguments. In the present paper, we 
suggest another, more general and justified condition, based on optimization procedure. 

3 Optimization method 

The idea of the method is to choose Ai = Aso that this choice would not influence much 
all other parameters Ai and rii. In other words, the factor approximant (3) should weakly 
depend on the chosen A. For small x, approximant (3) tends to 

fl{x) ~ fo{x) (x ^ 0) , (6) 

which does not depend on A, hence 

For large x, approximant (3) behaves as 

fl{x) r:^ fo{x)Bx^ (a; ^00), (7) 

where 

Nk Nk 

B = J{AT, /3 = E^- (8) 

i=l i=l 
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and Nk—{k + l)/2. Differentiating form (7) over A yields 

A, (11^ Bin. If). 

The solutions to Eqs. (5) are expressed through the unknown A = Ai, so that Ai = Ai{A) 
and Tii = ni{A). From Eqs. (8) it follows that B = B{A) and j3 = l3{A). Inverting 
these equations gives A — A{B) and A — Ai(l3). Equating A{B) = Ai{l3) implies that 
P = P{B), and, hence, /3 = P{B{A)). Therefore we have 

|j/,.w./„w.^(i + B.„.g^. (9) 

In order that the factor approximant (7) be weakly sensitive with respect to A, it sufficient 
to set 

dB , , 

This is the optimization condition defining A, which we shall employ. 



4 Illustration of method 

Many realistic problems are so complicated that only a few terms of series (2) can be de- 
rived by means of perturbation theory. Therefore, illustrating the suggested optimization 
method, we shall concentrate our attention on the case of short series (2), containing just 
a few terms. We shall find the optimized factor approximant 

f3'\x) = fo{x){l + Axni + A^xr , (11) 

in which A is prescribed by the optimization condition (10). And we shall compare form 
(11) with the nonoptimized variants of the factor approximants 

f;{x)=fo{x){l + xr{l + A,xr , (12) 

and 

f:{x) = foix){i+A,xr{i+A2xr' , (i3) 

in which all parameters are defined by the accuracy-through-order procedure, that is, are 
given by Eqs. (5). 



4.1 Logarithmic function 

To show that the method of factor approximants can approximate transcendental func- 
tions, let us start with a logarithmic function 

/(x) = - ln{l + x). (14) 
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Its expansion in series (2) yields 



fo{x) = 1 . 



n + 1 

The first several coefficients explicitly are 

1 1 1 1 

«i — ~ 2' ~ 3' ^3 — — 0,4 — -. 

For the optimized factor approximant (11), we find 

A = 0.65513, A2 = 0.94080, rii = -0.272, 713 = -0.440. 

Its large-variable behavior is 

f^P'ix) ~ 1.405x-°-^^' . 
This can be compared with the behavior of the nonoptimized approximant 

f;{x) ~ 1.260x-°-^6^ . 

Note that the direct consideration of the asymptotic series (2) for large x has no sense, 
leading to either —00 or +00, asx— s>oo. 

The factor approximants extrapolate the asymptotic series (2) to the region of large 
variables. As an example, let us consider the accuracy of approximants (11), (12), and 
(13) for X = 100. For x = 100, we have 

/3^*(100) = 0.052 , /3(100) = 0.058 , /4(100) = 0.054 . 

Comparing this with the exact value /(lOO) = 0.046, we find the corresponding percentage 
errors 

e (/s"^*) = 13% , e (/*) = 26% , s (/;) = 17% . 
As is seen, the optimized approximant is more accurate than the nonoptimized ones. 

4.2 Factorial function 

Factor approximants can well extrapolate even rather complicated transcendental func- 
tions, such as the factorial function 

1 ,vv/-/.r(i + i). (15) 



where r{x) is the gamma function. In expansion (2) for this function, one has fo{x) — 1 
and 

1 1 139 _ 571 



12' 288' 51840' 2488320 

The large-variable behavior of function (15) is 

f{x) ~ 0.398942a;°-^ {x ^ 00) . (16) 
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For the optimized approximant (11), we have 

A = 0.31614, = 0.31630, m = 474.945, ria = -474.418. 

The large- variable form of approximant (16) is 

f^P\x) ~ 0.419x°-^2^ (x^oo). 

This can be compared with the large-variable behavior of the nonoptimized approximant 

f;{x) ~ 0.119x°-^^^ (x^oo) . 

Again, the optimized factor approximant (11) essentially more accurately extrapolates the 
asymptotic series (2). For instance, ai x — 100, the error of the optimized approximant 
is about 15%, while the errors of the nonoptimized ones are several times larger. 

4.3 Debye-Hiickel function 

In the Debye-Hiickel theory of strong electrolytes [23] , there arises the function 

D{x) = ^ - ^ (l - ■ (17) 
The related series (2) contain fo{x) — 1 and the coefficients 

11 11 

^1 — ~ 77) 02 — 777, 03 — — — , 04 — 



3' 12' " 60' 360 

At large x, one has 

D{x) ~ 2x"^ (x^oo) . (18) 
For the optimized approximant (11), we get 

A = 0.2965, A2 = 0.2950, rii = 98.602, 712 = -100.221. 

Approximants (11), (12), and (13) extrapolate reasonably well series (2) up to x ~ 10. 
Comparing the accuracy of the approximants at x — 10, we get the percentage errors 

e (of) = -13% , e (d;) = -18% , e {DD = -5% . 

But for larger x ^ 10, the accuracy of all these approximants becomes not sufficient. 

4.4 Statistical integrals 

In the problems of statistical physics [23] , there occur the integrals of the following type 

POD g^" 

I(x) = / du . (19) 

^ ^ Jo 1 + xu ^ ^ 

Expanding this in powers of x, one has series (2), with fo{x) — 1 and the coefficients 
Oi = —1, 02 = 2, 03 = —6, a4 = 24. 



Such a series is strongly divergent. In the hmit of large x, the behavior of integral (19) is 
as {\nx)/x, similar to that of the logarithmic function (14). 
For the optimized approximant (11) we find 

A = 2.936042, A2 = 65.541043, rii = -0.340, ris = -1.559 x 10"^ 

To compare the accuracy of the extrapolation, let us consider the value of a; = 10. The 
corresponding percentage errors are 

e (if) = 56% , e (/*) = 1% , e {ID = 20% . 

In this case, the most accurate extrapolation is achieved by means of the factor approxi- 
mant Il{x). 

A close type of the integral, met in some problems of statistical physics [23] , reads as 

Its expansion in series (2) gives /o(a;) = 1 and the coefficients 

oi — 2, a2 — —2, 03 = —4, 04 = 24. 

This series are also strongly divergent. 

For the optimized approximant (11), we have 

yl = 3.370958 - 3.406509 i, A2 = A*, m = 0.119 + 0.175 i, n2 = n\. 

Comparing different approximants at the point x = 5, we find their percentage errors 

e (iT) = 12% , s {i;) = 5% , s (/:) = 18% . 

Here again, the most accurate extrapolation is provided by the approximant I^ix), al- 
though all three approximants (11), (12), and (13) have errors of the same order of 
magnitude. 

4.5 Partition function 

The typical structure of partition functions in statistical physics [23] is well described by 
the integral 

Z{g) = — / exp (-(^2 - gip') dip , (21) 
yTT J-00 ^ ' 

corresponding to the zero-dimensional field model, where g > ^ plays the role of the 
coupling parameter. Expanding Eq. (21) in powers of g yields series (2), with g instead 
of X, with /o(5') = 1, and the coefficients 

(-1)"^/ 1\ 
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The first coefficients expficitly are 



105 3465 675675 



^1 — ~ T) ^2 — T^TT) Q3 — 77^) ^4 — 



4' 32 ' 128 ' 2048 

The series is strongly divergent for any finite g. In the hmit of large g, one has 

Z{g) ~ 1.023^-°-25 (5 ^ 00) . (22) 

For the optimized approximant (11), we find 

A = 3.678882, A2 = 16.099756, m = -0.132943, ria = -0.016206. 

In the strong-coupling limit, we get 

Z"/\g) ^ 0.804^-°-^^^ (^ ^ 00) . 

In this limit, the approximants (12) and (13), respectively, are 

Z;{g) ^ O.mg-'-''' , Zl{g) 0.806g-' ''' (g ^ 00) . 

The accuracy of the approximants can be characterized by considering how well they 
approximate the amplitude and power of the strong- coupling limit, as compared to the 
exact behavior (22). We find that the amphtude is approximated with the percentage 
errors 

e (Z"^*) = -21% , e (Z;) = -10% , e (Z*) = -21% , 
while the power is predicted within the percentage errors 

s (^Zf) = -40% , e (Z*) = 38% , e {ZD = -48% . 

As is seen, the errors are of the same order of magnitude. Though these errors are not 
very small, we should not forget that the strong-coupling limit is predicted here on the 
basis of only weak-coupling terms, where g ^ 0, with no assumption on the large g ^ (yo 
behavior. Taking into account such a scarce information on just four weak-coupling terms, 
the prediction of the strong-coupling limit is not so bad. 

4.6 Anharmonic oscillator 

A great variety of problems in atomic and molecular physics, as well as in statistical physics 
and field theory, are modeled by the anharmonic oscillator, in which the anharmonicity 
corresponds to particle interactions. The anharmonic oscillator is represented by the 
Hamiltonian 

where x e (—00,00) and g e [0, 00). This model has been extensively studied by per- 
turbation theory with respect to the coupling parameter g. A detailed review article, 
describing the results of perturbation theory, is Ref. [24]. However, perturbation theory 
for the energy levels yields strongly divergent series for any finite value of the coupling 



9 



parameter g. Semiclassical approximation [24] can be used only for highly excited levels, 
but is not correct for low energy levels. In order to find effective limits of divergent series, 
the optimized perturbation theory [19,20], has been advanced, where the optimization is 
provided by introducing control functions governing the series convergence. The models 
of anharmonic crystals have been treated in this way [25-27]. Different variants of the op- 
timized perturbation theory for anharmonic oscillators have been considered [25-33] . The 
use of the optimized perturbation theory requires the introduction of control functions 
into the starting zero-order approximation. These control functions are to be defined at 
each approximation order, which makes calculations at high orders rather complicated. 
Our aim in the present paper is to show that the self-similar factor approximants can also 
be applied to the anharmonic-oscillator problem. The advantage of the factor approxi- 
mants is that they can be constructed on the basis of the pure pertubative series, without 
introducing control functions, which makes their use simpler. 

Calculating the ground state energy e{g) by means of perturbation theory [11,24] yields 
series (2) in powers of with 

k{g) = «o = 2 

and with the coefficients 

ai = ^, 02 = -2.625, 03 = 20.8125, 04 = -241.2890625. 

The strong-coupling limit of the ground-state energy is 

e{g) ~ 0.6679865/^/3 ^ 00) . (24) 

For the optimized approximant (11), we obtain 

A = 3.79054, A2 = 14.657819, Ui = 0.224, 712 = 0.044. 

In the strong-coupling limit, we have 

efig) ^ 0.759g'''' {g ^ 00) . 

Other approximants, in this limit, give 

el{g) ^ 0.611^'^-^^° , elig) c 0.755g' ''' {g ^ 00) . 

The amplitude of the strong-coupling limit is predicted by different approximants within 
the errors 

e (ef ) = 13% , e (e*) = -8.5% , e (el) = 13% , 
and the power of this limit is predicted within the errors 

e (ef ) = -19% , s (el) = 77% , e (el) = -31% . 

Again, we should not forget that these predictions for large g ^ 00 are based on the very 
scarce information of only four terms of the weak-coupling expansion, when g ^ 0. Going 
to the higher-order approximants essentially improves the accuracy. But let us recall that 
in the present paper we are interested in the possibility of extrapolating the weak-variable 
expansions to the large- variable range, when just a few first terms of such expansions are 
available. 
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4.7 Critical indices 



In this section, we demonstrate the apphcabihty of the method to calculating the critical 
indices related to the second-order phase transitions. We consider the </7^ - field theory, 
for which the critical indices can be calculated by means of the Wilson e-expansion, where 
e = 4 — d, with d being the space dimensionality (see, e.g., [34]). Let us be interested in 
the critical index u = i'{e), which, under the ^-expansion, is a function of s. One usually 
obtains the expansion for the function 

Treating e as an asymptotically small parameter, one gets [34] the expansion 

f{s) ~ 1 + ai£ + 026^ + a^e^ + a^e^ , (26) 
in which £ ^ and the coefficients are 



oi = -0.1665, 02 = -0.05865, 03 = 0.06225, 04 = -0.1535. 

However, the physical value of the index v corresponds to the space dimensionality c? = 3, 
hence, to the value s = 1. This means that it is necessary to extrapolate the asymptotic 
series (26) to the finite values of e. 

Accomplishing the extrapolation using the method of self-similar factor approximants, 
we extrapolate, first, the asymptotic series (26) getting f^ie)- Then we set e — 1. This 
allows us to define the corresponding approximation for the index 

^ (27) 



' 2/,*(l) • 

Following the scheme of Sec. 3, for the optimized approximant (11), we find 

A = 0.423, A2 = 0.425, ni = -306.868, 712 = 305.275. 

In different approximations, we obtain the values 

z/f * = 0.614 , 1/3* = 0.617 , = 0.634 . 

These values are very close to those derived by the Borel summation, resulting in 

u = 0.628 ± 0.001 [35] , u = 0.629 ± 0.003 [36] , 

and lattice numerical calculations 

1/ = 0.631 ± 0.002 [37] . 

Contrary to the latter numerical techniques, the method of self-similar factor approxi- 
mants is much more simple, giving the critical index of the close value of magnitude. 
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4.8 Polymer chain 

The expansion factor for a three-dimensional polymer chain with excluded-volume in- 
teractions is denoted as a{z), where 2; is a dimensionless coupling parameter [38,39]. 
Perturbation theory [38] results in an asymptotic series of form (2), with fo{z) — 1 and 
the coefficients 

ai = ^, a2 = -2.075385396, ag = 6.296879676, 04 = -25.05725072. 

Numerical fitting [38,39] shows that the polymer expansion factor can be represented by 
the empirical dependence 

a{z) = (1 + 7.524^ + 11.06^^)°'^^^^ (28) 

and in the strong-coupling limit, it behaves as 

aiz) ~ 1.531^°-=''^^ + 0.184z-°-5^^^ . (29) 

For the optimized approximant (11), we have 

A = 11.778451, A2 = 4.118898, rii = 4.84 x lO'^, = 0.31. 

The strong-coupling limit for this approximant is 

af{z) ~ 1.569^°-3i5 (^ ^ 00) . 

Approximants (12) and (13) give the strong-coupling limits 

al{z) ~ 1.424^°-^^° , al{z) ~ 1.560^°-3^° (^ ^ 00) . 

The strong-coupling limit here is of special importance, since it defines the strong- 
coupling exponent 

1 1 , Inafz) 

1/=- + - lim — 30 

2 4 2^00 1112; 

which is a kind of the critical index for the polymer chain [40] . It shows how the polymer 
expansion factor behaves at large coupling parameters, 

a{z) oc ^2(2-1) . 

For our approximants (11), (12), and (13), we find the index 

= 0.579 , u; = 0.610 , ul = 0.585 . 

This should be compared with the numerical estimates 

f = 0.5886 [38, 39], u ^ 0.5877 ± 0.0006 [40] . 

As is seen, the self-similar factor approximants yield the values that are very close to 
these numerical results. 
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5 Conclusion 



In the development of the method of self-similar factor approximants, we have suggested 
a solution to the problem of constructing the approximants of odd orders. We propose 
to define one of the parameters of the factor approximant by the optimization condition 
of form (10). The advanced technique is illustrated by a number of different problems 
with the mathematical structure typical of statistical physics and chemical physics. In 
these problems, it is necessary to extrapolate the asymptotic series, derived for an asymp- 
totically small variable, to the finite values of the latter and often one is interested in 
extrapolating the series to infinitely large variables. Calculations demonstrate that the 
method of self-similar factor approximants gives rather accurate results even when the 
number of terms in the asymptotic series is small. 

It is worth emphasizing that the method of self-similar factor approximants allows 
for accurate approximations not only for rational functions, as in the method of Pade 
approximants, but also for transcendental functions, as are considered in Sees. 4.1, 4.2, 
and 4.3. 

Moreover, there exists a wide class of functions that can be reconstructed exactly. This 
class of functions is formed by the following structures. Let Pn{x) be a polynomial in a 
real variable x of degrees n over the field of real numbers. And let ccj and (3j be arbitrary 
complex numbers. Compose the real-valued products of powers of such polynomials as 



n , n Pi 



The real values for these products are guaranteed, provided the complex powers ctj and 
Pj enter the corresponding products in complex conjugate pairs. Introduce the function 
Rmn{x) given by the ratio 



Rmn{x) 

where 



L\.i ^ n 



i j 

Denote the class of functions with the above structure through 

7^ = {Rmn{x) : M >0, N >0} . 

Suppose that a function f{x) is represented by an asymptotic series (2). Define the 
function 

= AW • 

Then the following theorem holds. 

Theorem. The function f{x), with expansion (2), is reproduced by means of the self- 
similar factor approximants (3), for any order k > M + N, exactly if and only if the 
function (p{x) pertains to the class TZ of functions of the form Rmn{x). 

The proof of this theorem is given in Rcf. [14] and also discussed in Ref. [15]. Concrete 
examples of the exact reconstructions are illustrated in Refs. [14-16,18]. 
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In conclusion, we may note that if the perturbative series do not have the form of an 
expansion in integer powers, as in Eq. (2), it is often easy to reduce them to that form 
by a change of variables. For instance, if we have an expansion in powers of y^, by the 
change x = g"^ we come to the expansion in integer powers for g. 
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